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ABSTRACT 

If MACHOs are black holes of mass ~ 0.5Mq, they must have been formed in the 
early universe when the temperature was ~1 GeV. We estimate that in this case in our 
galaxy's halo out to ~ 50kpc there exist ~5x 10 8 black hole binaries whose coalescence 
times are comparable to the age of the universe, so that the coalescence rate will be 
~ 5 x 10~ 2 events/year/galaxy. This suggests that we can expect a few events/year 
within 15Mpc. The gravitational waves from such coalescing black hole MACHOs can 
be detected by the first generation of interferometers in the LIGO/VIRGO/TAMA/GEO 
network. Therefore, the existence of black hole MACHOs can be tested within the next 
five years by gravitational waves. ( submitted to Apj Letters April 11 1997) 

Subject headings: gravitation — black holes — dark matter - — gravitational lensing - 
— Galaxy: halo 
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1. Introduction 

The analysis of the first 2.1 years of photometry 
of 8.5 million stars in the Large Magellanic Cloud 



(LMC) by the MACHO collaboration flAlcock et al 



1996a ) suggests that 0.62+q?, of the halo consists 



of MACHOs of mass 0.5±oJMq in the standard 
spherical flat rotation halo model. The preliminary 
analysis of four years of data suggests the existence 
of at least four additional microlensing events with 
tdur ~ 90days in the direction of the LMC ( Pratt 
1997 ). The estimated masses of these MACHOs are 
just the mass of red dwarfs. However, the contri- 



bution of the halo red dwarfs to MACHO events 



should be small since the observed density of halo 
red dwarfs is too low (Bahcall et al. 1994, Graff & 
Freese 1996a|, |Graff fc Freese 1996h| ). As for white 



dwarf MACHOs, the mass fraction of white dwarfs 
in the halo should be less than 10% since, assuming 
the Salpeter initial mass function (IMF), the bright 
progenitors of more white dwarfs than this would be 
in conflict with the number counts of distant galaxies 
( jCharlot fc Silk 1995| ). If the IMF has a sharp peak 
around 2Mfl, then the fra ction could be 50% or so 
( Adams fc Laughlin 1996] ), sufficient to explain the 
MACHO observations. The existence of such a pop- 
ulation of halo white dwarfs may or may not be con- 



sistent with the observed luminosity function (Gould 



1997, Lidman 1997, Freese 1997). In any case, future 



observations of hig h velocity wh ite dwarfs in our so- 
lar neighborhood ( Lidman 1997 ) will prove whether 
white dwarf MACHOs can exist or not. 

If the number of high velocity white dwarfs turns 
out to be large enough to explain the MACHOs, then 
stellar formation theory should explain why the IMF 
is sharply peaked at ~ 2Mq. If it is not, there arises 
a real possibility that MACHOs are absolutely new 
objects such as black holes of mass ~ 0.5Mq which 
could only be formed in the early universe, or bo- 
son stars with the mass of the boson ~ 10~ 10 eV. Of 
course it is still possible that an overdense clump of 



MACHOs exists toward the LMC flNakamura, Kan 
Y'a fc Nishi 1996|), MACHO s are brown dwarfs in the 



rotating halo flSpiro 1997|) or MACHOs are stars in 
the thick disk ( |Turncr 1997| ). 

In this Letter we consider the case of black hole 
MACHOs (BHMACHOs). In this case, there must 
be a huge number (at least ~4x 10 11 ) of black holes 
in the halo and it is natural to expect that some of 
them are binaries. In 52 we estimate the fraction 



f(a,e)dade of all BHMACHOs that are in binaries 
with semimajor axis a in range da and eccentricity 
e in de. We then use this distribution to estimate 
two observable event rates: First (end of §2) the rate 
of microlensing events we should expect toward the 
LMC due to binaries with separation ;> 2 x 10 14 cm; 
our result is in accord with the observation of one such 
event thus far (Bennett et al. 1996). Second (§3), 



the rate of coalescence of BHMACHO binaries out to 
15Mpc distance. The gravitational waves from such 
coalescences should be detectable by the first interfer- 
eomters in the LIGO/VIRGO/TAMA/GEO n etwork 



(Parish 1996| , [Brillct 1996| , [Tsubono 1996| , |Hough 
1996| ), and our estimated event rate is a few events 
§4 we discuss some implications of our 



per year. In 
estimates. 



2. Formation of Solar Mass Black Hole MA- 
CHO Binaries 

Since it is impossible to make a black hole of mass 
~ 0.5Mq as a product of stellar evolution, we must 
consider the formation of solar mass black holes in 



the very early universe (Yokoyama 1995, ledamzik 
1997). Our viewpoint here, however, is not to study 
detailed formation mechanisms, but to estimate the 
binary distribution that results. 

The density parameter of BHMACHOs, VIbhm, 
must be comparable to fib (or Vic dm) to explain the 
number of observed MACHO events. For simplicity, 
we assume that BHMACHOs dominate the matter 
energy density, i.e., fl — f2 bhm j though it is possible 
to consider other dark matter components in addition 
to BHMACHOs. To determine the mean separation 
of the BHMACHOs, it is convenient to consider it 
at the time of matter-radiation equality, t — t eq . At 
this time, the energy densities of radiation and BH- 
MACHOs are approximately equal and are given by 
p eq = 1.4 x 10~ 15 (flh 2 ) g/cm 3 , where h is the Hub- 
ble parameter in units of lOOkm/s/Mpc. Correspond- 
ingly, the mean separation of BHMACHOs with mass 
Mbh at this time is given by 

x = {M BH /p eq )z 

= 1.1 x 10 16 {M bh /MqY (Qh 2 )-hm. (1) 

We set the scale factor R to unity at t = t eq , so x can 
also be regarded as the comoving mean separation. 
Note that the Hubble horizon scale at t — t eq is L eq ~ 
1.1 x 10 21 (Vlh 2 ) ~ 2 cm. 



2 



During the radiation dominated era, the total en- 
ergy inside the horizon increases as R 2 . Since the 
Jeans mass in this era is essentially the same as the 
horizon mass, black holes are formed only at the time 
when the horizon scale is equal to the Schwarzschild 
radius of a BHMACHO. Thus the scale factor at the 
formation epoch becomes 



Rj = 



= 1.2X10-' 

C L e q 



(M, 



BH 



(2) 



The age of the universe and the temperature Tj at 
R = Rf are ~ 10 -5 sec and ~ 1 Gev, respectively. 

As a foundation for computing the distribution 
function /(a, e) for BHMACHO binaries, we assume 
that the BHMACHOs are created with a distribution 
of comoving separations x that is uniform over the 
range from an initial physical separation equal to the 
black hole size (which turns out to be so small that for 
the computations that follow we can approximate it 
as zero) to a maximum separation x — x. We also as- 
sume that the BHMACHOs' motions, if any, relative 
to the primordial gas have been redshifted to neglible 
speeds by the time their mutual gravitational attrac- 
tions become important. 

Consider a pair of black holes with mass Mbh and 
a comoving separation x < x. These holes' masses 
produce a mean energy density over a sphere the size 
of their separation given by pbh = peq{x / xR) 3 . This 
becomes larger than the radiation energy density p r = 

Peq/R* for 

R > R m = {x/xf . (3) 

This means that the binary decouples from the cos- 
mic expansion and becomes a bound system when 
R = Rm- Note that the background universe is still 
radiation dominated at this stage. 

If the motion of the two black holes is not dis- 
turbed, then the binary system cannot obtain any 
angular momentum so they coalesce to a single black 
hole on the free fall time scale. However, the tidal 
force from neighboring black holes gives the binary 
enough angular momentum to keep the holes from col- 
liding with each other unless x is exceptionally small. 

We refer to the semimajor axis and the semiminor 
axis of the binary as a and b, respectively, and we 
estimate a as 



d — X t^iix — X I X 7 



(4) 



We denote by y the comoving separation of the near- 
est neighboring black hole from the center of mass 



of the binary. Then b can be estimated as (tidal 
force) x (free fall time) 2 



b = 



GMbH xR m {xR m f 



{yRmf GM BH \ V/ 
Hence, the binary's eccentricity e is given by 



1 



(5) 



(6) 



Since (by assumption) x and y have uniform prob- 
ability distributions in the range x < y < 5, the prob- 
ability distribution of a and e is 



/(a, e) da de 



18x y x dx dy, 
(3/2)o*£-*e(l - e 2 yi dade(7) 



From the condition that y < x, the maximum value 
of the eccentricity for a fixed a is given by e max — 

J I — (a/x)i. Integrating /(a, e) with respect to e, 
we obtain the following distribution of the semimajor 
axis 



f a (a)da = - 



(!)'-©' 



- (8) 



From Eq. (g), it is found that the fraction of BH- 
MACHOs that are in binaries with a ~ 2 x 10 14 cm 
is - 8% and - 0.9% for Q,h 2 = 1 and 0.1, respec- 
tively. This estimated fraction of ~10 AU size BH- 
MACHO binaries is slightly smaller than the observed 
rate of binary MACHO events (one binary event in 8 
observed MACHOs), but the agreement is good given 
the statistics of small numbers. 

3. Gravitational Waves from Coalescing BH- 
MACHO Binaries 

We consider here short period BHMACHO bina- 
ries. Their coalescence times due to the emission of 



gravitational waves are approximately given by (Pe- 



ters and Mathews 1963, Peters 1964) 



-] (l-e 2 )*, 



where to — 10 year and 



a = 2 x 10 



li 



(MbhX 



(9) 



(10) 



is the semimajor axis of a binary with circular orbit 
which coalesces in to. Integrating Eq. (m) for a fixed 
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t with the aid of Eq. (^), we obtain the probability 
distribution for the coalescence time ft(t) as 



Mt)dt = - 



t 



dt , s 



where t max — to (x/oq) . If the halo of our galaxy 
consists of BHMACHOs of mass - 0.5M©, ~ 10 12 
BHMACHOs exist out to the LMC. The number 
of coalescing binary BHMACHOs with t ~ to then 
becomes ~ 5 x 10 8 for Qh 2 = 0.1 so that the 
event rate of coalescing binaries becomes ~ 5 x 
10~ 2 events/year/galaxy. If, however, the BHMA- 
CHOs extend up to half way to M31, the num- 
ber of coalescing binary BHMACHOs with t ~ to 
can be ~ 3 x 10 9 and the event rate becomes ~ 
0.3 events/year/galaxy. Both of these estimates are 
much larger than the best estimate of the event 
rate of coalescing neutron stars based on the statis- 
tics of binary pulsar searches in our galaxy, ~ 1 x 
10~ 5 events/year/galaxy (Phimiey 1991 , Narayan, Pi 



ran fc Shcmi 1991 , Van den Hcuvel fc Lorimcr 1996 ) 



The detectability of these waves by interferometers 
is most easily discussed in terms of the waves' "char- 
acteristic amplitude" h c (Eq. (46b) of Thorne 1987, 
with a well-known factor 2 correction): 



M °- 21 (^Sff '( 



1/6 / 



.WOHzJ \20Mpc 
(12) 

KereMchirp = (MiM 2 ) 3 / 5 /(Mi+M 2 ) 1/5 is the "chirp 
mass" of the binary whose components have indi- 
vidual masses Mi and M 2 . This h c is to be com- 
pared with an interferometer's "sensitivity to bursts" 
hsB = 11 [fSh(f )} 1 ^ 2 , where S/j(/) is the spectral den- 
sity of the interferometer's strain noise (cf. Eq. (Ill) 
of Thorne 1987, where hsB is denoted h 3 / yr ). This 
hsB(f) = hs/ yr (f) is plotted in various publications, 
e.g. Abramovici et. al. (1989) and Thorne (1995). It 
has a minimum (optimal sensitivity) at a frequency 
/ ~ 100 Hz. For the first LIGO and VIRGO in- 
terferometers, which are expected to be operational 
in 2001, that minimum is hsBmin — 3 x 10~ 21 . The 
GEO600 and TAMA interferometers, with their some- 
what shorter armlcngths, will have hsBmin & little 
worse than this in 2001. LIGO/VIRGO should be 
able to detect coalescing binaries, with high confi- 
dence, out to the distance for which h c = hsBmin 
at the optimal frequency / ~ 100 Hz. Inserting 
Mi = Mi ~ O.5M for BHMACHOs into the above 



equations, we see that the first LIGO/VIRGO inter- 
ferometers in 2001 should be able to see BHMACHO 
coalescences out to about 15Mpc distance, i.e., out 
to the VIRGO cluster, where our estimates (~ 1/100 
years in each galaxy like our own) suggest an event 
rate of several per year. 

LIGO R&D for the first interferometers is now 
nearing completion and is beginning to be redirected 
toward interferometer enhancements, for which the 
sensitivity goal is a factor 10 improvement, to hsBmin - 
3 x 10~ 22 (Barish et. al. 1996). In the mid 2000s, with 
these enhancements in place, LIGO should be able to 
see BHMACHO coalescences out to about 150Mpc, 
which would give a few events per year even if the 
event rate is 1000 times smaller than our estimates, 
~ 10~ 5 events/year/galaxy. 

4. Discussion 

In this Letter we have estimated the distribution 
function of binary BHMACHOs in order of magni- 
tude. It is possible to compute the distribution func- 
tion more accurately by N-body numerical simula- 
tions. This is an important, challenging numerical 
problem. 

Our estimated event rate for coalescing BHMA- 
CHO binaries is comparable to or greater than the 
most optimistic upper limit for binary neutron star 



coalescences (Phinney 1991), which are one of the 
most important sources of gravitational waves. Co- 
alescing neutron stars are also regarded as possible 



sources of the gamma-ray bursts (Mcszaros 1995). If 
so, then the detection of gravitational waves should be 
accompanied by a g amma-ray burst . If we consider 
the fire ball model ( Meszaros 1995 ) the time delay 
between the gravitational waves and the gamma-rays 
should be ~ 1 sec. By contrast, in the coalescence of 
binary BHMACHOs the emission of gamma-rays is 
not expected. This may enable us to distinguish co- 
alescing binary BHMACHOs from coalescing binary 
neutron stars. 

If gamma-rays are not emitted by coalescing binary 
neutron stars, we may still use their observed chirp 
masses M c hi rp to distinguish them from BHMACHO 
binaries. The chirp mass can be measured from the 
gravitational waves to a fraction of a per cent accu- 
racy, which is much less than the expected spreads 
of BHMACHO masses and neutron star masses. The 
masses of neutron stars in binaires are expected to be 
~ 1.4Mq corresponding to a chirp mass of 1.2Mq. 
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This is supported observationally as well as theoreti- 
cally since the mass of the iron core before the collapse 
is ~ IAMq. 

Therefore if the chirp masses of BHMACHO bi- 
naries are much smaller than IMq it is possible to 
identify them. However if the mass of a BHMACHO 
is ~ IMq, coalescing BHMACHO binaries and coa- 
lescing neutron star binaries may be indistinguishable 
from the detected gravitational waves. In principle, 
however, black holes absorb some of the gravitational 
waves so that the evolution of the binary in the last 



few minutes (Cutler et al. 1993) is slightly different 
from that of a neutron star binary of the same mass. 



Th : difference arises in the 2.5th post-Newtonian or- 



der for the Kerr black hole case (Tagoshi 1997 ) and 
4t h order for the Sch warzschild black hole case ( |Pois- 
son fc Sasaki 1995 ). This difference might be de- 



tectable, though it will present a challenging problem 
for gravitational wave data analysis. 

It is known that the Silk damping of density pertur- 
bations on small scales causes distortions of the CMB 
spectrum by dumping acoustic energy into heat and 
thenc e into the CMB ([Daly 199l| , |B arrow fc Colc£ 
19S0, |Hu, Scott fc Silk 1994|) . Thus it is important to 
examine if the large primordial density perturbations 
that are needed for BHMACHO formation are com- 
patible with the observed upper limit of CMB spectral 



distortions ( [Mather et al. 1994| ). 

For a rms amplitude of density perturbations 5(1) 
on a comoving scale I in the radiation dominated 
era, the fraction of the universe that turns into black 
holes is given by f(l) w 5(1) exp (-1/I85(l) 2 ) (|Carr 
197|5]). In the present case, f(l) is equal to Rf given 
by Eq. (||) and I is equal to the comoving scale Ibh 
of the density perturbations that give rise to BHMA- 
CHOs. Recalling that we fixed the scale factor R 
to unity at the time t eq of matter-radiation equality, 
the present comoving scale corresponding to Ibh is 
given by (1 + z eq )l BH = (1 

5x 10 17 (Mbh/Mq) ' 

at t = t eq . Since Rf ~ 10 -8 , we need 5(Ibh) 2 ~ 
4 x KT 3 . 

On the other hand, Hu, Scott and Silk (1994) 
showed that the observational limit constrains the rms 
amplitude of density perturbations at horizon entry 
to be S(l) 2 < 10~ 4 for I > Id, where Id is the comov- 
ing scale of the Silk damping at the double Comp- 
ton thermalization time, which at present is given by 



Z e q)Rj 1 GM B HC 



cm, where z eq is the redshift 



terest Ibh is approximately 3 orders of magnitude 
smaller than the scale Id to which the observational 
constraint applies. 

Correspondingly, if we assume the primordial den- 
sity perturbation spectrum to have a power-law shape, 
P(k) oc fc™, on scales greater than Ibh, the above re- 
sult implies that n J> 1.6. Interestingly, this value 
is consistent with the one suggested by the observed 
CMB anisotropics (Bennett et al. 1994) and it is 
marginally consistent with the COBE normalization 
([Hu, Scott fc Silk 1994]). It may be that this rather 



blue spectrum can be produced by a variant of the so- 



called hybrid model of inflation(Garcia-Bcllido, Lindc 



Wands 1996) 
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